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1 0 Formation of normal control subject data base 
20 Patient Selection and screening 

101 Patient Presents to a Clinic with Memory 
Complaints, Family History of Alzheimer's 
Disease, ApoE4 Positive or Related Reasons 

i 

102 Patient Presents to a Clinic with History of 
Substance Abuse, Concern over Possible Brain 
Damage or Related Reasons 

i 

103 Patient Presents to a Clinic with a History of Head 
Injury, Cognitive Complaints, Headaches, Blurry 
Vision or Related Reasons 

104 Patient Presents to a Clinic with Family History of 
Parkinson's Disease, Movement Disorder or 
Related Reasons 



110 Standard Clinical Evaluation Performed: Baseline 
Laboratory Tests; Computerized Tomography or 
Magnetic Resonance Imaging Performed to 
Investigate Potential Causes of Complaint; and 
Neuro-psychological Testing is Performed 

115 Other Causes of Cognitive Decline Ruled Out 
(Neoplasm, Endocrine Imbalance, Infection, 
Nutritional Deficiency etc. 

V 

120 Patient is Referred to the PET Center for Evaluation 
of Cognitive Decline Index 

y 

FIG. 1 Normal Subject Controls and Patient Selection for Development of the 

Cognitive Decline Index 
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200 
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225 



Patients are injected with 10 mCi of FDG through a peripheral 
intravenous line. 

V 

Control subjects were injected with 5 mCi. 

V 

Patient/Subject Scanned; Raw Brain Image Data Collected 
in 2D Mode, Septa in Place 

V 

10 Minute Emission Scan is Obtained for Patient 

V 

5 Minute Transmission Scan is Obtained for Patient 

V 

20 Minute Emission Scan is Obtained for all Controls 

V 

5 Minute Transmission Scan is Obtained for all Controls 

y 

Ordered Subject Expectation Maximization (OSEM) 
Reconstruction 

V 

Segmented Transmission Scan Attenuation Correction 

Y 

Addition of Z-axis filter 



FIG, 2a Derivation of the Processed Digital Brain Image 
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230 Images are then reformatted for clinical evaluation on the 
GE Advance Workstation (Sun™ Ultra 60), and a Copy 
Prepared for Export to the Research Workstation for Use 
with SPM 



235 Copy Prepared for Export to the Research Workstation for 
Use with SPM 



240 Scans Converted to Analyze 7.5™ (Analyze Direct, 
LENEXA KS) Format 



t 

245 Images Are then further Processed and Analyzed with 
SPM99 (SPM, Friston et al, 1995a) and Implemented in 
Matlab 



V 



250 Management of the Image Data is Achieved by Utilizing 
the Digital Imaging and Communication in Medicine 
(DICOM) Format 



■V 

260 Convert to Signed, 16-bit, "Byte-swapped" Image 



265 "Byte-swapping" Option Taken Depending on Workstation 
Type 



270 Image Transformed into Talairach Space 



280 Image Spatially Transformed using Gaussian Kernel 



FIG. 2a (Cont'd) Derivation of the Processed Digital Brain Image 
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300 SPM Group Analysis: Patients Compared to Normal 
Controls 

310 Extrema in Changes in Metabolic Brain Activity Evaluated 



320 SPM(Z) Maps Created; Tables of Regional Maxima 
Compiled 



330 Convert to Normal Distribution Z-score 

340 Identification and Location oi Clusters with Highest 
Significance Using Significance Thresholds Uncorrected 
for Multiple Comparisons ^ 

350 Coordinates of Locations of Significant Interest Identified 
from SPM results listing all maxima greater than 8mm apart 

360 Specific Loci from this Analysis are Used as Centers for 3-D, 1 
cm Diameter Spherical (VOIs) 

370 Identification of Specific Loci in Brain as Volumes of 
Interest (VOI) created with the Marsbar® plug-in (Bret et al, 
2002) for SPM 



380 intensity of Each of the Voxels within the Spherical VOI is Read 
and the Average is Obtained 

390 Weights (Set GDI,) for each VOI are Calculated Based on the 
Frequency of Abnormality of the VOI Data from All the Study 
Patients as Compared to All Controls with Greater Weight Given 
to Increased Frequency as Determined by Intensity Range 
Overlap for each Volume of Interest Ratio 

i 

FIG. 2b Derivation of Region Location and Identification of VOIs 
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391 Import VOI Data Into Spreadsheet 

i 

392 Determine Intensity Range Overlap for each VOI Ratio 

393 Create Weights for each Intensity Extreme 

i 

394 Create Weighted VOI Ratio 

395 Scale and Normalize Ratio 



400 Calculation of GDI: 



CDI^C x +{ If fx tin ]/[ YjV&lm ] 

Where Xj denotes the j th Increased Intensity Value; 

Vj denotes the j* Weight for the j* Increased 
Intensity Value; 

Yi denotes the 1 th Decreased Intensity Value; and 

Wi denotes the i* Weight for the i* Decreased 
Intensity Value. 

C x is the correction factor used to normalize the 

dataset. 



i 



410 Weights of Set CDli are then used as a baseline for calculation 
of a second set of Weights (Set CDI 2 ) to calculate CDI 2 . Set 
CDI2 is calculated by iterative optimization of each weight to 
maximally separate the patient from the controls 



FIG. 2b (Cont'd) Derivation of Region Location and Identification of VOIs 
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500 GDI is Compared to Established Normal Range of Values; 
Presence of Normality or Abnormality is Determined 

J 

510 If GDI Reading is Negative, Patient Educated About the 
Clinical Course of Potential Illnesses and Signs to Watch 
out for; Potential Benefit of Preventative Measures such as 
Antioxidants, Brain Exercises, Beneficial Diet and 
Adequate Rest are Discussed 

I 

520 If GDI Reading is Positive, Patient Educated about the 
Meaning of a Positive GDI Reading, Projected Clinical 
Course of Illness, Benefit of Medication, Potential Benefit 
of Preventative Measures such as Antioxidants, Brain 
Exercises, Beneficial Diet and Adequate Rest are Discussed 

I 

530 Following either Step 620 or Step 630, Results are Given to 
Referring Physician and Patient Scheduled for Re- 
evaluation in one year. 

if 

540 Patient Data is Stored in the Comprehensive Patient 
Database 



FIG. 2c Patient Diagnosis and Clinical Recommendations 
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Fig. 2D Construction of the Neural Network Based GDI (CDInn) 

600 Append the Neural Network Toolbox to the Matlab Path 

I 

610 Type the command "nntool" at the Matlab command prompt to start the NN graphical 
user interface 

620 Configure an artificial neural network with the features and parameters: 
Name: CDInn 

Type: Feed-forward back-propagation 

Input ranges: [-1 1; -1 1; -1 1; -1 1; -1 1; -1 1; -1 1; -1 1; -1 1;] 

Training function: trainlm 

Adaption learning function: Learndgm 

Performance function: MSE 

Layers: 2 

Number of neurons: 2 

Transfer function neuron 1: purelin 

Transfer function neuron 2: tansig 

630 VOI regions obtained from the specified locations: 
Regions of increased activity: 
vermis, 
motor, 
R pons, and 
cerebellar nuclei 

Regions of decreased activity: 
posterior cingulate, 
L parietal, 
R parietal, 
L tempo ralOl, and 
L tempora!02 

4- 

640 Intensity data from the above regions are entered into the NN in the order: 
vermis, 
motor, 
R pons, 

cerebellar nuclei, 
posterior cingulate, 
L parietal, 
L tempo ralOl, 
L temporaI02, 
R parietaL 

650 VOI Data are mean-normalized by subject, i.e. divided by the average of the 9 regions 
measured for each subject. This gave each subjects data set a mean value =1. 

4, 

660 The datasets are next divided into comparably sized subgroups (cl denotes control 
subgroup # 1, pi = patient subgroup #1, etc): [cl=16 ? c2=17, pl=16, p2=16]. 

4, 
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670 Four mixed groups were made from these by pairing controls with patients to give: [dpi, 
clp2, c2pl,c2p2]. 

680 These datasets are then exported (from Excel) to a text file, subsequently imported into 
Matlab and converted to a .mat file for use by Matlab. 



690 Rescaling is performed to ensure the input range is appropriately matched for the tansig 
function on the output layer of the NN, which will classify the input data in the range ± 1 . 
This is carried out in the main Matlab window at the command prompt. As an example, 
the dataset clpl is rescaled using the premnmx function: 



700 The rescaled output from this function is labeled and saved as clpl_rescale.mat. This is 
done for each of the four paired datasets listed above. 

710 These four files are then used as "inputs" for neural net training. For training of the 
neural net, each input file had a "target" file to which the neural network tried to match 
trie input dataset to. These files were labeled to match the dataset from which they were 
used, thus for clplrescale.mat, the "target" file was labeled clplt.mat. Each target file 
consisted of a single row of digits with the same number of columns as the input file. For 
each patient (column) in the clpl_rescale dataset, a "-1" was present in the clplt.mat 
file, and for each control, a "1" was assigned. Thus, the four input and target datasets 
were: i 



[clpl_rescale,minp,maxp]==premnmx(clpl) 




Input Target 

clpl_rescale.mat clplt.mat 
clp2_rescale.mat clp2t.mat 
c2pl_rescale.mat c2plt.mat 
c2p2_rescale.mat c2p2t.mat 



720 



730 



8/17 



WO 2005/094686 



PCT/US2005/010439 



Inputs Weights 
i 



Bias Transfer Output 
b function 




Nei%vork 
inputs 



Hidden 



Output Network 
layer output 




F&.3A 



Z30 



9/17 



WO 2005/094686 



PCT/US2005/010439 




Fig .3B 
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Design matrix 



set-levet 


cluster-level 






voxeH 




ve\ 




x 1 y,z{mrn} 


P o 










T 








0.173 14 


0.000 


1945 


O.OOO 


0.106 


5.29 


< 


4.47) 


0.000 


-4 


-70 30 










0.4B3 


4.54 


i 


3 .97) 


0.000 


-14 


-68 16 










0.616 


4.37 


< 


3.86) 


O.OOO 


-4 


-53 28 




o.ooo 


2610 


o.ooo 


0.788 


4.16 




3 .70) 


0.000 


-42 


-74 36 










0.604 


4.13 


< 


3 .68) 


0.000 


-56 


-56 16 










0.681 


4.01 


( 


3 .39) 


0 .000 


-60 


-56 -8 




0.862 


202 


O.064 


0.683 


4.01 


{ 


3.59) 


0.000 


-6 


14 -24 




0.247 


450 


0.009 


0.951 


3.84 


( 


3.47) 


0.000 


52 


-64 36 










1.000 


2 .96 


( 


2.7?) 


0.O03 


54 


-50 44 




1.000 


sa 


O.3 02 


0.998 


3 .49 


( 


3.20) 


0.001 


-48 


IS -22 




0.991 


114 


Q.153 


0.998 


3 .48 


( 


3.19) 


0.001 


-4 


-12 10 










1.000 


2.81 


< 


2 .64) 


0.004 


-8 


-30 4 




0.996 


101 


0.177 


0.999 


3.43 


( 


3.15) 


0.001 


4 


38t -16 




0.944 


162 


O.093 


1. 000 


3.37 


( 


3 .10) 


0.001 


-10 


14 -2 










1.000 


2 .90 


( 


2 .71) 


0.003 


-12 


0 12 




0.9&1 


131 


Q.127 


1.000 


3.32 


( 


3.07) 


0.001 


-20 


-i4 -2a 










1.000 


2.91 


i 


2 .72) 


0.003 


-30 


-22 -20 




0.8&8 


191 


0.071 


1.000 


3.22 


( 


2.98) 


0.001 


-38 


20 -2 




0.997 


95 


0.189 


1.000 


3.21 


( 


2.97) 


0.001 


-26 


48 16 

.36 IB 




0.993 


110 


0.160 


1.000 


3.05 


( 


2. 84) 


0.002 


-2 










1.000 


2 .69 


< 


2 .54) 


0.006 


-2 


42 2 




1.000 


69 


0.261 


1.000 


3.00 


( 


2.80) 


0.003 


-24 


-38 -6 




1.000 


SO 


0.338 


1.000 


2.93 
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2.74) 


0.003 


-40 


46 6 



^ table shows at most local maxima > 8.0mm apart per cluster 

HelipTT^ Degrees of freedom = [1 .0. 33.01 



Extirif threshold: k = 50 voxeis,~p~= &%38~(1 .000 corrected) Smoothness FWHM == 13.5 13.7 16,0 {mm) = 6.7 6,9 8.0 Jvoxels} 
Expected voxels per cluster, <k> = 58.849 Search .volume: S = 1815544 rrirn A 3 = 228943 voxels = 559.8 resell 



Expected number of clusters, <c> = 10.48 



Voxel size: [2.0, 2.0, 2.0] mm (1 resel = 370.95 voxels) 



Statistics: 



smgie clus^ summ corrected for entire volume) 



■iflnrTflfifort 



voxel-level 
~T 



■.Mnfinmftfitari 



x I y 1 z{mrn} 



o.ooo 



2610 



o.ooo 



no. v 



0.788 
0.804 
0.SB1 
0-892 
0.305 
0. 93 S 
0.977 
0.987 
0.993 
1.000 
1.O00 
1.000 
1.000 
1.000 



4.16 
4.13 
4.01 
3 .99 
3 .96 
3.89 
3 .74 
3 .86 
3. SO 
3 .03 
3 .01 
2.99 
2 .83 
2.76 



( 3.70) 

< 3.60) 
{ 3.59) 

< 3.58) 
( 3.56) 

< 3.50) 

< 3.39) 
{ 3.33) 
( 3.21) 

< 2.83) 



< 2.B1) 

< 2.79) 

< 2.66) 
{ 2.60) 



0.000 
O.OOO 
0.000 
0.000 
O.OOO 
0.000 
O.OOO 
0.000 
0.O01 
0.O02 
0.O02 
0.O03 
0.O04 
0.005 



-42 
-56 
-60 
-62 
-64 
-58 
-50 
-42 
-52 
-56 
-54 
-62 
-62 
-32 



-74 3 6 

-56 16 

-56 -8 

-36 -6 

-30 -22 

-46 -22 

-60 42 

-62 44 

-62 28 

-46 36 

-68 16 

-24 -8 

-18 -24 

-56 40 
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Increases in MCI 
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Statistics: 

"se^eveT*'" 1 ' 



volume summary (p-va!ue& corrected for entire volume 



contrasts) 





Design matrix 



rv 

4e 



cluster-level 



voxei 
"7 



ievei 



nnnmriRfitRTf 



x,y,z{rnrTt} 



0.003 



0.002 



0.313 
0.003 



0.001 
0.530 
0.6 S3 



745 



188 
678 



873 
138 
110 



0.000 



0.026 
0.000 



0.000 
0.051 
0.077 



0.027 
0.528 
0.773 
0.120 
0.202 
0.622 
0. 925 
0.379 
0.434 
0.536 
0.337 
0. 937 
0. 994 
0.995 
1.000 



5.86 
4 .48 
4 .18 
5.23 
5.00 
4.37 

3 .92 
4.66 

4 .60 
4 .40 
4. OS 
3 .89 
3 .59 
3.57 
3 .30 



( 4.82) 

< 3.93) 

< 3.71) 
( 4.43) 

< 4.28) 
{ 3.85> 

< 3.52) 

< 4.07) 

< 4.02) 

< 3.87) 
( 3.65) 
( 3.50) 
( 3.28) 
( 3.26) 
( 3.05} 



0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 

o.oao 

0.000 
0.000 
o.ooo 

0.001 

0.001 

0.O01 



-16 -24 52 
-26 -26 66 
-6 -34 62 
26 -54 64 
10 -22 -3 0 
14 -38 -34 
-6 -28 -24 
34 -16 68 
26 -32 52 
14 -38 68 
62 12 36 
62 O 22 
62 O 34 
-6 -54 -16 
4 —SO -6 



table shows at most focal maxima > B.Omm apart per cluster 



Height threshold: T = 3.00, p = 0,003 (1 .000 corrected) Degrees of freedom = [1 XX 33.0] 

Extent threshold: k = 100 voxels, p = 0.090 (0.739 corrected) Smoothness FWHM = 13.5 13.7 16.0 {mm! = 6.7 6.9 8.0 -{voxels} 
— - • - Search volurne:S = 1815544 mm A 3 = 226943 voxels = 56b,8 resels 



Expected voxels per cluster, <k> = 35.683 
Expected number of clusters, <c> = 1 .34 



Voxel size: [2.0, 2.0, 2.0] mm (1 resel = 370.95 voxels) 



Fig. 5 
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Fig. 8 
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